Observation of the Cabibbo-suppressed decay — ► pK n + 
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We report the first observation of the Cabibbo-suppressed charm baryon decay H^" — > pK~ff ■ 
We observe 150 ± 22 events for the signal. The data were accumulated using the SELEX spec- 
trometer during the 1996-1997 fixed target run at Fermilab, chiefly from a 600 GeV/c S~ beam. 
The branching fractions of the decay relative to the Cabibbo-favored — > T, + K~tv + and 
H+ -> H-7T+7T+ are measured to be B(E+ -> pK-ir+)/B(E+ -> E+ K~ tt+) = 0.22 ± 0.06 ± 0.03 
and S(H+ -> pK~-K + )/B(E+ -> H _ 7r + 7r+) = 0.20 ± 0.04 ± 0.02, respectively. 
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The study of Cabibbo-suppressed (CS) charm decays 
can provide useful insights into the weak interaction 
mechanism for non-leptonic decays [jlj. The observed 
final state may arise either from direct quark emission at 
the decay stage or, in some cases, from quark rearrange- 
ment due to final-state scattering. By comparing the 
strengths of CS decays to their Cabibbo-favored (CF) 
analogs, one can, in a systematic way, assess the contri- 
butions of the various mechanisms. 

Modern methods for calculating non-leptonic decay 
rates of the charm hadrons employ heavy quark effective 
theory and the factorization approximation . Nonethe- 
less, the three-body decays of charm baryons are pro- 
hibitively difficult to calculate due to the complexity of 
associated final state interactions. Measurements of the 
relative branching fractions of charm baryon states, both 
CF and CS, give additional information about the struc- 
ture of the decay amplitude and the validity of the fac- 
torization approximation. 

Many CS hadronic decays of charm mesons have been 
measured by both e + e~ collider and fixed target exper- 
iments. Until now the only CS charm baryon decay 
reported with significant statistics is A+ — ► pK~K + 
and its resonant state, A+ — » pcf) |||]]. This transition 
requires internal W-emission and may be inhibited by 
color-suppression effects. In this paper, we present the 
first observation of the CS decay S+ — > pK~ir + and de- 
termine the branching fractions of this decay relative to 
the CF S+ -> T.+K-TT+ and S+ -► S-7T+7T+ modes. All 
three of these modes involve external W-emission. 

Figure [l] shows a spectator diagram with external W- 
emission for S+ decaying into a CF and a CS mode. (The 
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FIG. 1. An example of spectator diagrams for ; 
(a) CF H+ -> and (b) CS H+ -> pK~ V 



decays, 



other CF S~ mode interchanges s and u quark lines and 
produces a dd pair from the vacuum instead of a uu pair.) 
The decay process is similar in the two modes except for 
the flavor change in the weak decay Cabibbo-Kobayashi- 
Maskawa matrix element (V cs vs. V c d). Based on a sim- 
plified dimensional analysis, we would expect -B(S+ — > 



pK-n+)/B{E+ -> Y.+K-K+) = 2.1 x a x tan 2 6» c where 
8 C is the Cabibbo angle and a, of order unity, measures 
the ratio of the hadronic contributions to the weak decay 
matrix element. The factor 2.1 reflects the relative phase 
space of the two modes. By comparing the a factors for 
the entire family of CS non-leptonic charm baryon de- 
cays one may develop a picture of the dominant features 
of the amplitudes. 

The SELEX(E781) experiment at Fermilab is partic- 
ularly well-suited for studying CS charm baryon decays 
because of its excellent particle identification capabilities. 
SELEX is a high energy hadroproduction experiment us- 
ing a 3-stage spectrometer designed for high acceptance 
for forward interactions [xp — 2p^ AI /y/s > 0.1). The 
experiment emphasizes understanding charm production 
in the forward hemisphere and the study of charm baryon 
decays. Using both a negative beam (~ 50% E~, ~ 50% 
7T-) and positive beam (~ 92% p, ~ 8% tt+), SELEX 
recorded 15.2 billion interaction events during the 1996- 
1997 fixed target run. The majority particles (£~, tt~) in 
the 600 GeV/c negative beam and the protons in the 540 
GeV/c positive beam were tagged by a beam transition 
radiation detector. The data were accumulated from a 
segmented target (5 foils: 2 Cu, 3 C, each separated by 
1.5 cm) whose total thickness was 5% of an interaction 
length for protons. 

The spectrometer had silicon strip detectors to mea- 
sure the beam and outgoing tracks, giving precision pri- 
mary and secondary vertex reconstruction. Transverse 
position resolution (a) was 4 /im for the 600 GeV/c beam 
tracks. The average longitudinal vertex position resolu- 
tion was 270 fim for primary and 560 /im for secondaries, 
respectively. Momenta of particles deflected by the an- 
alyzing magnets were measured by a system of propor- 
tional wire chambers (PWCs), drift chambers and sili- 
con strip detectors. Momentum resolution for a typical 
100 GeV/c track was o p jp w 0.5%. The absolute mo- 
mentum scale was calibrated using the Kg mass. For 
D° — ► K~tt + decays the average mass resolution was 
9 MeV/c 2 independent of D° momenta from 100-450 
GeV/c. Charged particle identification was done with a 
Ring Imaging Cerenkov detector (RICH) ||, which dis- 
tinguished K ± from 7r ± up to 165 GeV/c. The proton 
identification efficiency was greater than 95% above pro- 
ton threshold (sa 90 GeV/c). For pions the total mis- 
identification probability due to all sources of confusion 
was less than 4%. 

Interactions were selected by a scintillator trigger. The 
trigger for charm required at least 4 charged tracks after 
the targets as indicated by an interaction counter and at 
least 2 hits in a scintillator hodoscope after the second 
analyzing magnet. It accepted about 1/3 of all inelastic 
interactions. Triggered events were further tested in an 
on-line computational filter based on downstream track- 
ing and particle identification information. The on-line 
filter selected events that had evidence of a secondary 
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vertex from tracks completely reconstructed using the 
forward PWC spectrometer and the vertex silicon. This 
filter reduced the data size by a factor of nearly 8 at 
a cost of about a factor of 2 in charm written to tape 
(as determined from a study of unfiltered Kg, A and 
D° — ► K~ir + + c.c. decays). Most of the charm loss 
came from selection cuts that were independent of charm 
species or kinematic variables and which improved the 
signal/noise in the final sample. 

Results presented here come from the first processing 
pass through all data. In this analysis, secondary vertex 
reconstruction was attempted when the x 2 P er degree of 
freedom (x 2 /dof) for the fit of the ensemble of tracks to 
a single primary vertex exceeded 5. All combinations of 
tracks were formed for secondary vertices with as many 
as 5 prongs. Secondary vertices were tested against a 
reconstruction table that specified selection criteria for 
each charm decay mode. Secondary vertices which oc- 
curred inside the volume of a target or after the first 
plane of the silicon vertex detector were rejected. 

The pK~TT + sample was selected by the following re- 
quirements: (i) primary and secondary vertex fits each 
have x 2 /dof < 5; (ii) vertex separation significance 
L/cr > 8, where L is the longitudinal separation between 
primary and secondary vertices and a is the error on L; 
(iii) the reconstructed momentum vector from the sec- 
ondary vertex points back to the primary vertex with 
X 2 /dof < 4; (iv) the momentum of the tt + > 5 GeV/c; (v) 
£(K)/£(tt) > 1 for K identification and C(p) /£(tt) > 1 
for p identification, where C is the mass-selection like- 
lihood function for the RICH; and (vi) the sum of the 
squared transverse momenta of the three daughter par- 
ticles with respect to the parent direction, Ep 2 ^ > 0.3 
(GeV/c) 2 _ 

Figure |] shows the first observation of the CS H+ — > 
pK~ir + decay mode. The inset of the figure shows the 
invariant mass distribution from the entire selected mass 
range of reconstructed pK~Tr + candidates. The large 
peak in the inset is the decay A+ — * pK~ir + . The bump 
at the right is the CS S+ — > pK~ir + decay. The dis- 
tribution has an artificial high-mass cutoff because of a 
maximum mass cut imposed in this first-pass analysis. 
Monte Carlo simulations show that the observed num- 
ber of signal events is systematically reduced only by 
(2 ± 3)% due to this cutoff. The background character- 
istics are very similar both in the number per mass bin 
and lifetime through the entire mass region, so that the 
background subtraction under the S+ peak is stable. The 
lifetime both of background and signal events for the 
region has been investigated. The background apparent 
lifetime is quite short, while the S+ lifetime is consis- 
tent with the Particle Data Group average ||. A de- 
tailed lifetime analysis is underway and will be reported 
elsewhere. The mass is determined using a Gaussian of 
variable centroid, width and area. The number of signal 
events is determined by a sideband subtraction method. 
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FIG. 2. The invariant mass distribution of pK~n + . 



The signal window is determined from the A+ data and 
is set at 45 MeV/c 2 . The background is a second order 
polynomial fit to the entire mass region excluding a 65 
MeV/c 2 wide window centered at each mass peak. The 
5+ yield is the difference between the summed events in 
the 45 MeV/c 2 signal window centered at 2.467 GeV/c 2 
and the extrapolated background. Variations in the ex- 
clusion window size change the extrapolated background 
by fewer than 2 events. There are 150 ± 22 events for 
S+ -> pK-ir+ at a S+ mass of 2467.4 ± 1.3 MeV/ c 2 . 
The statistical significance for the signal, S/ 'V S + B , is 
6.9 ± 0.6 in which S is the number of signal events and B 
is the number of background events in the signal region. 

Background contributions from charm baryon (S" and 
f2°) decay to the CS decay are expected to be negligi- 
ble; they produce an even number of secondary vertex 
tracks and their decay lengths are short compared to the 
minimum flight path cut. Excellent proton identifica- 
tion strongly suppresses background from charm meson 
decays (three-body D + and Df ). 

One normalizes the CS mode by comparing it to CF 
modes that have a hyperon in the final state. In the 
first reconstruction pass, hyperons (S ± ,S _ ) were iden- 
tified only inclusively in a limited decay interval (5-12 
m downstream from the target): candidates had p > 40 
GeV/c and no hits assigned along the trajectory in the 14 
chambers after the second analyzing magnet. This cate- 
gory of tracks gives unique S + identification but is am- 
biguous between E~ and 3~. Figure || shows two CF S+ 
modes decaying to S~7r + 7r + and S + K~n + , respectively. 
The shaded areas in Fig. 3(a) and 3(b) are the estimated 
reflection from A+ — > £~7t + 7t + and A+ — > £ + 7r~7r + , re- 
spectively. The shapes are determined by a Monte Carlo 
simulation and the areas are normalized to the observed 
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number of signal events in the A+ data. For this part 
of the analysis, additional cuts were applied to all data. 
The pion momentum threshold was raised from 5 to 10 
GeV/c and the transverse component of the momen- 
tum with respect to its line of flight was required to be 
less than 0.3 GeV/c. 
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FIG. 3. Charm baryon E^ signals with hyperon partial re- 
construction, (a) E~7r + 7r + and (b) T, + K~tt + . The solid line 
is the background fit. 

The total acceptance (geometrical acceptance and re- 
construction efficiency) for decay modes of interest was 
estimated by embedding Monte Carlo charm decay tracks 
into data events. Momentum and energy were not con- 
served in the process, but studies indicate this has little 
effect on the single-charm acceptance calculation. Events 
were generated with an average transverse momentum 
(pt) = 1-0 GeV/c and longitudinal momentum distribu- 
tions as observed for the data. Detector hits, including 
resolution and multiple Coulomb scattering smearing ef- 
fects, produced by these embedded tracks were folded 
into the hit banks of the underlying data event. The new 
ensemble of hits was passed through the SELEX off-line 
software. The acceptance is the ratio of the number of 
events in the signal to the number of embedded events in 
a particular mode. 

The method was checked by measuring two well- 
determined relative branching fractions, B(D° — > 
K'tt' 
K-K+ir+)/B(D 



K--K + Tz + ir-)/B{D° 



B(D^ 



^ = 1.94 ± 0.07±g;?i and 



K--K+TT+) = 0.093 ± 
0.010±g 

^006 wnere the first error is statistical and the sec- 
ond is systematic. Systematic errors include the differ- 
ence of the ratio for charge-conjugate states, and uncer- 
tainties of signal yields and acceptance estimations. The 
results agree well with the world average Q . For the D + 
decays, the vertex significance cut (L/a) was increased to 





Accept ance(%) 


Events 


Ej — > E 7T + 7T + 


2.46 ± 0.04 


130 ± 15 


E+ - 


1.09 ±0.02 


53 ±10 


E+ -» pK~n + 


7.08 ± 0.06 


76 ±13 



TABLE I. Summary of observed events and estimated ac- 
ceptance for modes using the tighter cuts. The errors are 
statistical only. 



20 to suppress background from Df — > K~ K + ir + . This 
costs little signal, since the lifetime of D + is 2.4 times 
longer than that of . 

Using this procedure, we measure the relative branch- 
ing fraction for S+ — > Y, + K~ir + and S+ — > S~7r + 7r + to 
be (0.92 ± 0.20 ± 0.07). The systematic error is due to 
uncertainties in the background subtraction and relative 
acceptance estimation. This result is comparable to the 
CLEO measurement, (1.18±0.26±0.17) @. The number 
of events and estimated acceptance for the three observed 
E!+ modes with the common set of cuts are summarized 
in table |. 

The branching fraction of the CS decay — » pK~n + 
relative to the CF S+ — * Y, + K~ir + is measured to be 
0.22 ± 0.06 ± 0.03 which corresponds to (4.2 ± 1 .1 ± 0.5) x 
tan 2 c . The a factor is 2.0±0.5±0.2. To put this m some 
perspective, the corresponding a factor for the CS decay 
A+ — * pK~K + relative to the CF mode A+ — > pK~n + 
is 2.5 ± 0.6, using the world average for the branching 
fraction and correcting for phase space. Thus, the fac- 
tors for these two different baryons are similar. In con- 
trast, the CS ratios of D mesons depend strongly on fi- 
nal state multiplicity, suggestive of sizable final state in- 
teractions. For the E~tt + it + mode, B(S+ — > pK~ir + )/ 
B(S+ -> E-TT + n + ) is measured to be 0.20 ± 0.04 ± 0.02. 
Systematic errors include uncertainties in (i) the relative 
acceptance estimation, (ii) background subtractions due 
to the reflection from other hyperon modes, and (iii) the 
pK~Tr + yield due to the imposed mass cutoff. 

In summary, we observe the CS decay — > pK~n + 
at a mass of 2467.4 ± 1.3 MeV/c 2 with 150 ± 22 sig- 
nal events. The relative branching fractions of the 



the CF S+ -> v - 
measured to be 
7T+) = 0.22 ±0.06 ±0.03 and 
-► S-7T+7T+) = 0.20 ± 0.04 ± 



Y.+K- 



decay S+ — > pK tt + to 
and — > S~7r + 7r + are 
pK-ir+)/B(Z+ -> E+K 
B(S+ - pK-ir+)/B(E+ 
0.02, respectively. 

The authors are indebted to the staffs of Fermi Na- 
tional Accelerator Laboratory, the Max-Planck-Institut 
fur Kernphysik, Carnegie Mellon University, and Peters- 
burg Nuclear Physics Institute for invaluable technical 
support. This project was supported in part by Bun- 
desministerium fur Bildung, Wissenschaft, Forschung 
und Technologie, Consejo Nacional de Ciencia y Tec- 
nologfa (CONACyT), Conselho Nacional de Desenvolvi- 
mento Cientffico e Tecnologico, Fondo de Apoyo a la In- 



4 



vestigacion (UASLP), Fundagao de Amparo a Pesquisa 
do Estado de Sao Paulo (FAPESP), the Israel Science 
Foundation founded by the Israel Academy of Sciences 
and Humanities, Istituto Nazionale de Fisica Nucleare 
(INFN), the International Science Foundation (ISF), the 
National Science Foundation (Phy #9602178), NATO 
(grant CR6.941058-1360/94), the Russian Academy of 
Science, the Russian Ministry of Science and Technology, 
the Turkish Scientific and Technological Research Board 
(TUBiTAK), the U.S. Department of Energy (DOE 
grant DE-FG02-91ER40664 and DOE contract number 
DE-AC02-76CHO3000), and the U.S.-Israel Binational 
Science Foundation (BSF). We also acknowledge a useful 
conversation with Austin Napier. 



Present address: Imperial College, London SW7 2BZ, 
U.K. 

t Present address: Instituto de de Ffsica da Universidade 

Estadual de Campinas, UNICAMP, SP, Brazil 
* Deceased 

§ Present address: Dept. of Physics, Wayne State Univer- 
sity, Detroit, MI 48201 

** Present address: Universitat Freiburg, 79104 Freiburg, 
Germany 

ft Present address: Physik-Department, Technische Univer- 
sitat Miinchcn, 85748 Garching, Germany 
Present address: Instituto de Ffsica Teorica da Universi- 
dade Estadual Paulista, Sao Paulo, Brazil 
§§ Present address: Lucent Technologies, Naperville, IL 
*** Present address: Max-Planck-Institut fur Physik, 

Miinchen, Germany 
tt1 " Present address: Motorola Inc., Schaumburg, IL 
Mt Present address: Carnegie-Mellon University, Pittsburgh, 
Pittsburgh, PA 

§§§ p resen t address: Deutsche Bank AG, 65760 Eschborn, 
Germany 

[1] J.G. Korner and M. Kramer, Z. Phys. C55, 659 (1992). 
[2] M. Bauer, B. Stcch and M. Wirbel, Z. Phys. C34, 103 
(1987). 

[3] CLEO Collaboration, CP. Jessop et al, Phys. Rev. Lett. 

B423, 56 (1994). 
[4] E687 Collaboration, P.L. Frabetti et aJ., Phys. Lett. 

B314, 477 (1993). 
[5] J. Engelfried et al., Nucl. Instrum. Methods A 431, 53 

(1999). 

[6] Particle Data Group, C. Caso et al, Eur. Phys. J. C3, 1 
(1998). 

[7] CLEO Collaboration, T. Bergfeld et al., Phys. Lett. 
B365, 431 (1996). 



5 



